Experimental Observations of Group Synchrony in a System of Chaotic Optoelectronic Oscillators 
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We experimentally demonstrate group synchrony in a network of four nonlinear optoelectronic oscillators 
with time-delayed coupling. We divide the nodes into two groups of two each, by giving each group different 
parameters and by enabling only inter-group coupling. When coupled in this fashion, the two groups display 
different dynamics, with no isochronal synchrony between them, but the nodes in a single group are isochronally 
synchronized, even though there is no intra-group coupling. We compare experimental behavior with theoretical 
and numerical results. 



The last years have seen a vast increase in the interest in 
coupled dynamical systems, ranging from few coupled ele- 
ments to complex networks |Q][2]]. Besides the focus on net- 
work structure and topology, the research area of synchroniza- 
tion in networks has grown rapidly ||3] |4). The groundbreak- 
ing work on the master stability function (MSF) by Pecora and 
Carroll has bridged the gap between topology and dynamics 
by allowing predictions about synchronization based solely on 
the nodes' dynamics and the eigenvalue spectrum of the cou- 
pling matrix 0. 

While the MSF theory was originally developed for identi- 
cal, isochronous synchronization, more complex patterns of 
synchronization are observed in applications in, e.g., neu- 
ral systems, genetic regulation, or optical systems I6 HT31 . 
These patterns include, for example, sublattice synchroniza- 
tion in coupled loops of identical oscillators with heteroge- 
neous delays [ 16 1, pairwise synchronization of pairwise iden- 
tical nodes coupled through a common channel [ 17 1, and more 
general group synchronization 1 18 1. In group synchronization 
the local dynamics in synchronized clusters can be different 
from the dynamics in the other cluster(s), which extends the 
possibility of synchronization behavior to networks formed of 
heterogeneous dynamical systems, as they appear in a vari- 
ety of applications. Moreover, these synchronous patterns can 
be observed even when there is no intra-group coupling. Sor- 
rentino and Ott have generalized the MSF approach to group 
synchronization |18|, and recent work by Dahms et al. con- 
siders time-delayed coupling of an arbitrary number of groups 
ED. 



In this Letter, we demonstrate the successful realization 
of group synchronization of chaotic dynamics in an array 
of four optoelectronic oscillators. Optoelectronic oscillators 
with time-delayed feedback have been found to show a mul- 
titude of different dynamical behaviors ranging from steady- 
state to chaotic dynamics depending on parameters Il20lj25l . 
In this work we experimentally demonstrate group synchrony, 
where the two groups display different fluctuation amplitudes. 
Remarkably, the two groups of synchronized oscillators are 
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FIG. 1 . (Color online) (a) Schematic of four nodes separated into two 
groups, A (red, solid) and B (blue, dashed), (b) Experimental setup 
of a single node, showing coupling to the other nodes according to 
the configuration in (a). 



not directly coupled to each other; they are only coupled to 
those of the other groups. 

The experimental setup consists of four optoelectronic 
feedback loops, which act as the four nodes of the network. 
We consider several coupling schemes. In the first one, the 
nodes are coupled together in the configuration shown in 
Fig.[TJa) in order to form two groups. There are no direct cou- 
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pling links between two nodes in the same group. However, 
a node is coupled bidirectionally to both of the nodes in the 
other group. In this experiment, the coupling strength, e, and 
coupling delay, r, are the same for all coupling links. How- 
ever, the parameters of the nodes differ depending on which 
group the nodes are in. Both of the nodes in group A are iden- 
tical, and both of the nodes in group B are identical, but the 
nodes in group A are not identical to the nodes in group B. 
In Fig.[TJa), the coupling links are shown in black (arrows in 
each direction to indicate bidirectional coupling), and the self- 
feedback of the nodes is indicated by the gray (colored) lines 
and arrows. 

A schematic of a single node is shown in Fig. |TJb), where 
red lines indicate optical fibers, and black or green lines indi- 
cate electronic paths. In each node, light from a diode laser 
passes through a Mach-Zehnder modulator (MZM), whose 
output light intensity is cos 2 (a: + (f>o) for an input voltage sig- 
nal x. There is a controllable bias phase of the MZM, which 
we set to be 00 = f . The optical signal is split into three 
equal signals: one is the feedback signal, and the other two 
are the coupling to the two nodes in the opposite group. A 
photoreceiver converts the feedback optical signal to an elec- 
trical signal, which is one of the two inputs to the DSP (digital 
signal processing) board. The incoming optical signals from 
the two nodes of the other group are combined optically be- 
fore a second photoreceiver converts the composite coupled 
signal to an electronic signal, which is the second input of 
the DSP board. The DSP board implements the feedback and 
coupling time delays, which are the same for this experiment 
(t = 1.4 ms), and a diffusive coupling scheme. The feed- 
back signal is scaled by a factor of 1 — e, while each incoming 
signal to a node is scaled by a factor of ~, for the global cou- 
pling strength, e, and the number of links incoming to a node, 
n m . For the configuration shown in Fig.[TJa), n m = 2 for all 
nodes, but in general, n; n can be different for each node. The 
feedback and coupled signals are summed on the DSP board. 

The DSP board also implements a digital filter, which is 
a two-pole bandpass filter with cutoff frequencies at 100 Hz 
and 2.5 kHz and a sampling rate of 24 kSamples/s, and also 
scales the combined signal by a scale factor, which controls 
the feedback strength, which we denote /3. The output of the 
DSP board is amplified with a voltage amplifier, whose out- 
put drives the MZM. Although (3 is a combination of gains of 
the photoreceiver, amplifier, and other components, the DSP 
board is the only place where the gain is changed. 

For this experiment, all parameters except for (3 are identi- 
cal in all four nodes. We keep (3 identical among the members 
of each group but allow a different f3 for each group, denoted 
by @W and ft B l Pr evious studies have revealed the wide 
variety of behaviors that are possible for this type of system, 
depending on the value of (3 [22 1. For this study, we have used 
a range of (3 from to 10, with the experiments focusing on 
cases of f3 > 3, for which the system displays chaos (with 
some periodic windows) when the nodes are not coupled. 

For each run of the experiment, the nodes are started from 
random initial conditions. This system has a time delay, so the 
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FIG. 2. (Color online) Maximum Lyapunov exponent A max as a func- 
tion of /3' A ' and /3' B ': (a) in the longitudinal directions 7 = ±1, 
(b) in the transverse direction 7 = 0. White areas correspond to 
A max = 0. Black dot indicates values of and /?' s ' used in this 
experiment. 



initial condition will be a function of time. Thus, we record 
the random electrical activity at the input to the DSP in the 
absence of coupling and feedback for 1 second to provide the 
initial states for the nodes. After recording an initial condition, 
we enable feedback for 4 seconds, which is long enough for 
transients to disappear. At the end of this period, we enable 
coupling. Data are taken after transients have died out. 

The system of coupled feedback loops can be well- 
described by a mathematical model with a system of time de- 
lay differential equations for the voltages input to the MZMs 

x\ e 



and the vectors describing the states of the filters 
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where m and m! 7^ m denote the groups A or B, and i indi- 
cates the node within a group. E = ( ), F = 
( ), and G = ( 1 ) are constant matrices that describe the 
filter. The filter parameters are chosen as ujl — 2tt ■ 2.5 kHz 
and ljjj = 2-7T • 0.1 kHz. For a bipartite network with no intra- 
group coupling, we define the inter-group coupling matrices 
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where K is the overall coupling matrix for the entire net- 
work. For the configuration shown in Fig.[TJa), i,j = 1,2, 
and K( A ) = K' 5 ' = \ ( \ \ ) so that 
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Equations ([TJ and Q can describe the dynamics of the uncou- 
pled nodes if we set the coupling strength e = 0, as the second 
term in Eq. (|2]i represents the diffusive coupling scheme. 



3 



K 








K (fl > 






K 



K 



K 




time, / (ms) 



FIG. 3. (Color online) (a) Simulated motion in the synchronization 
manifold, obtained by numerically solving Eqs. l|6j and l|7j, show- 
ing the predicted group-synchronous state. (b),(c),(d) Experimen- 
tally measured time traces from three different network configura- 
tions (indicated by the coupling scheme and the coupling matrices 
K) that achieve group synchrony. All three networks have the same 
eigenspectra, but the configuration in (b) is symmetric while those in 
(c) and (d) are not. 
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FIG. 4. Correlation functions of 3.9 s of experimental data, (a) Au- 
tocorrelation functions for the dynamics of group A (left) and group 
B (right), with no coupling in the system, (b) Cross-correlation func- 
tions between the two nodes in group A (left) and group B (right), 
for coupled nodes, (c) Cross-correlation functions between one node 
in group A and one node in group B for the uncoupled system (left) 
and coupled system (right). 



Numerical simulations use a discrete-time implementation 
of these differential equations, as described in Ref. Il22l . The 
simulations of uncoupled and coupled systems are in excel- 
lent agreement with the experimental results for the variety of 
dynamical behaviors that can be observed. 

We will now investigate the existence and stability of the 
group synchronous solution, i.e., we will derive analytical 
conditions determining whether such a solution (in which the 
two nodes of each group are identically and isochronously 
synchronized, but there is no identical synchrony between 
nodes of different groups) exists for given values of and 
j3^ B \ and if it does, if that solution is stable. We use the ap- 
proach described in |[T8l[T9l . The condition for the existence 
of the group synchronous solution for a particular coupling 
configuration is that 

J2Kjf=c^,m = {A,B}, (5) 

3 

i.e., that the row sum of the matrices K*™ 1 ' is constant. For 



the work reported here, we fix &- A ' = = 1. 

The group synchronized dynamics for group m is given by 

ui m) (i) = Eu<™)(£) - cos 2 (4 m) (t - t) + <t> ), (6) 

x<THt) = G[(ui m \t) + e(ui m '\t) - u<T\t))]. (7) 

Linearizing Eqs. ([TJ and |2]i about the synchronous solution 
Us (to = A, B), we obtain the master stability equations: 

Su^(t) =E<5u M (<) - F/?( m > sin(2a;< m >(t - r) + 2<fo) 

x G[(l - e)Su {m \t -t)+ e 7 u (m,) (i - r)]. 

(8) 

In Eq. ([8), 7 is a parameter that is chosen from the eigen- 
value specttum of K. The largest Lyapunov exponent as a 
function of this parameter 7 is called the MSF. For the con- 
figurations presented here, the nonzero eigenvalues of K are 
1 and -1, and any remaining eigenvalues are zeros. There- 
fore, the stability results will be identical for any two-group 
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network whose nodes are described by Eqs. ([Tji and |2]) and 
whose coupling matrix is in the form of OJ, satisfies ([5j, and 
has identical rows for either K.^ or K^' (for a proof, see 
supplemental material). 

The eigenvalues 7 = — 1 and 7 = 1 in the master stability 
equation ([8} correspond to perturbations parallel to the syn- 
chronization manifold. The corresponding value of the MSF 
determines the dynamical behavior inside the synchronization 
manifold and is shown in Fig. |2|a) in dependence on the pa- 
rameters and ft B \ Negative, zero, and positive values 
denote fixed-point, periodic, and chaotic dynamics, respec- 
tively. Due to the inversion symmetry of the MSF for two- 
group synchronization [18, 19], the MSF values are identical 
for 7 = — 1 and 7=1. 

Transverse stability of the synchronization manifold is de- 
termined by using the eigenvalue 7 = in Eq. Q. Figure|2|b) 
shows the largest Lyapunov exponent in the transverse direc- 
tion, which is negative for almost the entire range of 13^ and 
that is shown, indicating that we expect the group syn- 
chronous solution to be stable for most parameters. 

To observe group synchrony in this system, we select dis- 
similar values of and (3^ B \ as shown by the black dots in 
Fig. [2] The global coupling strength is chosen as e = 0.8. The 
experimental values for p( A ' and were adjusted using the 
DSP board. The values of f}( A ) and used in simulation 
were established by varying the values close to the experimen- 
tal values to find nearby values which match best the dynami- 
cal behavior of the experiments for uncoupled nodes, obtained 
from the shape of the reconstructed attractor in phase space. 
Since the values determined experimentally as f}( A > — 7.6 and 
f3( B ) = 3.3 are subject to measurement uncertainties, it is not 
surprising that we find slightly different values in simulation, 
i.e., = 7.66 and f3^ B ' > = 3.28. Comparison of uncou- 
pled and coupled time traces in experiment and simulation is 



shown in the supplemental material, Fig. SI 



Figure [3] shows experimental and simulated time traces of 
the coupled system. The simulated traces in Fig. [3ja) show 
the behavior of any two-group system displaying stable group 
synchrony according to Eqs. |6} and (mi, with the parameters 
we have used here. Figure J3jb) shows experimental results 
for a system coupled according to Fig.[TJa). These time traces 
show that there is identical, isochronal synchrony between 
x± and a4 , and between x^ and x^\ but not identi- 
cal synchrony between the groups. Thus, this is an example 
of group synchrony. We also performed experiments on two 
asymmetric four- node configurations. These configurations 
were created by removing links from the original structure of 
Fig. [TJa), while preserving the constant row sum and eigen- 
values (1,-1,0, and 0) of K, keeping all other parameters the 
same. Their topologies and dynamics are shown in Figs.[3jc) 
and|3jd). Because these schemes are also described by Fig. [2] 
they also display group synchrony. In the experimental time 
traces, there are slight differences between the two traces of 
one group, due to the intrinsic experimental noise and mis- 
match we expect in any real system. An example of a larger 
network that displays the same behavior is presented in the 



supplemental material, Fig. S2 

To further examine the nature of the synchrony of this sys- 
tem, we calculate the correlation functions of the experimen- 
tal time traces, as shown in Fig. |4] for the topology shown in 
Fig.[TJa). For two variables y(t) and z(t), which each have a 
mean of zero, we define the correlation function C as a func- 
tion of time lag At 



C(At) 



(y(t)z(t + At)) 



(9) 



Figure |4|a) shows the autocorrelation functions for one node 
in each group when the nodes are uncoupled. The autocor- 

f A) 

relation of x\ shows only a peak at zero time lag, which 

( B) 

indicates chaotic dynamics, while the autocorrelation of x\ 
shows periodic dynamics, with correlation peaks at intervals 
of the time delay r = 1.4 ms. In Fig. |4|b), we show the 
cross-correlation functions of x^' with x^\ and of x[ B ^ 

( B) 

with x 2 for the coupled system, which confirms identical, 
isochronal chaotic synchronization between the two nodes in 
a single group. Figure |4jc) shows the cross-correlation func- 
tions between two nodes in different groups, without and with 
coupling. The uncoupled case has no correlation, as we ex- 
pect, but the coupled case has a high correlation peak at a 
lag of At = —1.4 ms. From this, we can see that there is 
time-lagged phase synchrony between the two groups, with 
the dynamics of group B leading the dynamics of group A 
by the system delay, r. However, the amplitudes of fluc- 
tuations of the two groups are still different after coupling, 
so there is no complete synchronization, and we have an in- 
teresting situation of the simultaneous coexistence of intra- 
group isochronal identical synchrony and time-lagged phase 
synchrony between the groups. 

In conclusion, we have examined a four-node system of 
nonlinear optoelectronic oscillators in the case where there 
are two groups of nodes with dissimilar parameters. Our ex- 
periments display the phenomenon of group synchronization, 
and we analyze the stability of the group synchronized solu- 
tions for chaotic dynamical states. It is remarkable that, al- 
though the coupling is entirely between the different groups 
and not within the groups, identical isochronal synchroniza- 
tion within each group is induced by this coupling, while the 
two groups are not mutually amplitude synchronized, as pre- 
dicted by our stability analysis using the generalized master 
stability function |[T8l[T9l . Thus the nodes of group B act as 
a kind of dynamical relay [27 1 for the nodes of group A, and 
vice versa. These results have been experimentally demon- 
strated with three coupling configurations, and conditions for 
observing group synchrony in other networks have been dis- 
cussed. 

Our observations go beyond previous work on sublattice 
and cluster synchrony, where the experiments focused on op- 
tical phase synchronization for coupled lasers without self- 
feedback J9j[l0j. Group synchronization in larger networks is 
a significant challenge for future experimental investigation. 
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SUPPLEMENTAL MATERIAL 



Figures [3jb),[3jc), and[3jd) show that stable group synchrony is experimentally observed for three different coupling configu- 
rations. Here we show that our stability analysis and the numerical computations in Fig. |2]apply to all of these coupling schemes 
and, more generally, to a whole class of networks, characterized by an arbitrary number of nodes in both the groups A and B. 

We define Na and Ng the number of nodes in group A and B, respectively. Then the couplings are fully described by the 
Na x Nb coupling matrix *KS A \ whose entries {K^} represent the intensity of the direct interaction from system j in group 

B to i in group A and the Nb x Na matrix T£( B \ whose entries {K^} represent the intensity of the direct interaction from 
system i in group A to j in group B. 

First we note that the motion in the synchronization manifold (Eqs. |6]) and |7|) applies to any network described by Eqs. ([T| 
and |2]), as long as the entries along the rows of the matrices and K^ B ) sum to one. Thus in what follows, we will limit our 
attention to the case that 



C (A) = C (B) = 



= 1. 



(SI) 



If assumption < |S 1 J > is verified, it follows that the maximum Lyapunov exponent of the synchronous solution shown in Fig. |2ja) 
does not depend on the details of the underlying network structure. 

According to Ref. 1 18 1, a master stability function approach to group synchronization is possible for any network described 
by Eqs. ([T} and |2]), under the assumption (SI I. For any such network stability depends on the eigenvalues of the matrix 



K = 



K(^) 
K( fi ) 



(S2) 



From Ref. 1 18 1 we see that these eigenvalues are 



A= [0,0,...,0](J[±VA 1 ,±VA 2 ,...,±v / Aiv mi J, (S3) 
where [0, 0, 0] denotes | Na — Nb \ zeros and A = {Ai, Ajv mi „ } denotes the spectrum of the N m i n x N m i n matrix, 



D 



K(^)K( S \ if Na < Ng, 



K( B )K( A \ if N B <N A , 
We proceed now under the assumption that either one of the two following conditions is satisfied, 

fa\ 



(S4) 



K (B) = 



W 

b 

w 



(S5a) 



(S5b) 



where a (b) is any Nb -dimensional row-vector (A^-dimensional row-vector) with its entries summing to one. Then the matrix 
D is in the form 



D = 



(d\ 

d 

\"d) 



(S6) 



where d is an N m i n -dimensional row -vector with sum of its entries equal one. Note that the underlying assumption is that either 
one of the matrices and K^ B ) is in the form of Eq. ( S5 1 (not necessarily both). Then the eigenvalue equation for the matrix 
D reduces to the following equation 

\Vi=d-V=£ i = 1, N min , (S7) 
where v = [v±,V2, —r v N min ] is an eigenvector, A is the associated eigenvalue, and i is a constant that does not depend on i. 



There are two possible solutions for Eq. ( S7 1 
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(i) I = 0. Then, A = and v is any vector orthogonal to d. There are (N m i n — 1) such vectors. 

(ii) i ^ 0. Then, v\ = V2 — ■■■ — VN min and from the sum of the entries of the vector d being equal one, it follows that A = 1. 
Hence by using Eq. ( |S3| l, and assuming satisfaction of either one of the Eqs. ( |S5a| i or ( |S5b[ ), the eigenvalues of the matrix K 

are 

a=[o,o,...,o]|J[-ML (S8) 

where here [0, 0, 0] denotes \Na + Nb — 2| zeros. It follows that stability of the group synchronous solution for any network 
that satisfies either Eqs. ( |S5[ a)) or ( |S5[ b)) is described by the plot in Fig. |2|b). 
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FIG. SI. Comparison of time traces from experiment (left column) and simulation (right column). Nodes in group A (B) are indicated by 
the red and solid (blue and dashed) lines, (a) Nodes are uncoupled, and the uncoupled nodes are completely unsynchronized. (b) Nodes are 
coupled according to the configuration in Fig. Ill a). There is identical, isochronal synchrony between x[ A ^ and x^ A \ and between x[ B ^ and 

( B) — 

x 2 , but not identical synchrony between the groups, so this is an example of group synchrony. In the simulation, the two traces in one 
group are exactly synchronized and are indistinguishable, as we expect from a simulation without noise or mismatch. In experiment, slight 
differences between synchronized traces arise from experimental noise and mismatch in the real experimental system. In both the experiment 
and simulation, the dynamics of the nodes in group B have a significantly smaller amplitude than those in group A, with qualitatively similar 
dynamics between simulation and experiment. 
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FIG. S2. Simulation of a seven node network displaying group synchrony. The coupling scheme and the associated coupling matrix are shown 
above the simulated time traces. The parameters are the same as those presented in the letter. Coupling is enabled at time t = 0, and the nodes 
become group synchronized so that all three red, solid lines of Group A are identically synchronized, and all four blue, dashed traces of Group 
B are identically synchronized. Note that the coupling matrix for the specific network configuration shown here satisfies Eq. ( |S5b[ >. Hence, 
the coupling matrix has eigenvalues -1, 1, and degenerate zeros. The stability of the group synchronous solution of this system is described by 
Fig.[2jb). The dynamics of the seven node, group synchronized system, are predicted by Fig.[2ja), as illustrated by Fig.|3[a). 



